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Context. Both the absolute mass-loss rates and the mechanisms that drive the mass loss of late-type supergiants are still not well 
known. Binaries such as a Sco provide the most detailed empirical information about the winds of these stars. 

Aims. The goal was to improve the binary technique for the determination of the mass-loss rate of a Sco A by including a realistic 
density distribution and velocity field from hydrodynamic and plasma simulations. 

Methods. We performed 3D hydrodynamic simulations of the circumstellar envelope of a Sco in combination with plasma simula- 
tions accounting for the heating, ionization, and excitation of the wind by the radiation of a Sco B. These simulations served as the 
basis for an examination of circumstellar absorption lines in the spectrum of a Sco B as well as of emission lines from the Antares 
nebula. 

Results. The present model of the extended envelope of a Sco reproduces some of the structures that were observed in the cir- 
cumstellar absorption lines in the spectrum of a Sco B. Our theoretical density and velocity distributions of the outflow deviate 
considerably from a spherically expanding model, which was used in previous studies. This results in a higher mass-loss rate of 
(2 ±0.5) x 10~ 6 M yr -1 . The hot Hn region around the secondary star induces an additional acceleration of the wind at large distances 
from the primary, which is seen in absorption lines of Ti n and Cr n at - 30 km s~ 1 . 
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1. Introduction 

In the late stage of the stellar evolution stars suffer a substan- 
tial mass loss by a continuous wind. This phenomenon plays 
a key role for the final development of individual stars and for 
the enrichment of the interstellar medium. However, the physi- 
cal mechanisms that govern the outflow of late-type giants and 
supergiants are not established. Possible driving mechanisms in- 
clude Alfven waves, a coustic waves, shock waves, and pulsa - 
tions (see the reviews of Laf on & Berruverll991l:IWillsonl2 000). 
As it is hard to distinguish between different driving mechanisms 
from a theoretical point of view, there is still a need for empirical 
studies to provide constraints. 

The most accurate mass-loss rates are obtained from obser- 
vations of eclipsing binaries of f Aur or VV Cep type. These 
systems consist of a late-type supergiant and a main-sequence 
companion that can be used as a "natural satellite" to probe 
the circumstellar e nvelope. For the b asic concepts of the binar y 
te chnique see e g. iReimersI dl987l) . lHack & Sticklandl (1 19871) . 
or lBaade) dl998l). Observa tions with hi gh spectral resolution of 
i Aur dBaade et al.lll996h and a Sco (Baa de & Refiners! l2007t 
Reime rs et al.l 2008) have shown that the circumstellar envelopes 
of these stars cannot be described by a continuous, spherically 
expanding model. One of the aims of this study was therefore to 
investigate the influence of the binary companion on the density 
and velocity distribution of the supergiant wind. 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, un- 
der program ID 076.D-0690(A), and on observations made with the 
NASA/ESA Hubble Space Telescope (program # 5952), obtained at the 
Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. 



The a Sco system (Ml lb + B2.5 V) is of particular impor- 
tance in this context, because its angular separation amounts to 
nearly 3", which allows th e observation of the e xtended enve- 
lope with spatial resolution. Rei mers et al.l d2008l) presented the 
most recent extensive observations, which were carried out with 
the Ultraviolet and Visual Echelle Spectrograph (UVES) on the 
Kueyen telescope of the Very Large Telescope array (VLT) and 
provide spatially resolved emission spectra of the Antares neb- 
ula. They showed that the mass-loss rate of a Sco A can be deter- 
mined by measuring the spatial extent of the Ha emission. They 
determined a mass-loss rate of M = (l.O5+O.3)xlO -6 M yr~' as 
a best match between the measured Ha distribution and plasma 
simulations, which were based on a spherically expanding model 
of the circumstellar envelope and led to a nearly spherical Hn 
region around the secondary that is rotationally symmetric with 
respect to the line connecting the two stars. They reported an 
apparent increase of the mass-loss rate towards larger distances 
from the supergiant. Moreover, they noted that the observed 
[Fe n] emission cannot be explained in the framework of their 
simplified model. 

Mid-infrared observations at the Keck n telescope of the dust 
in the circumstellar envelope of a Sco indicate a non-uniform 
dust distribution that ap pears to be rela ted to discrete ejections 
of mass from a Sco A dMarshet al.l200ll) . This is consistent with 
interferometric observati ons that were carri ed out at the William 
Herschel Telescope by Tuth ill et al.l d 19971) and revealed asym- 
metric structures ("hotspots") at the surface of the supergiant. 

Baade & Reimers (2007) presented spectra of a Sco B, 
which were obtained with the Goddard High-Resolution 
Spectrograph (GHRS) on the Hubble Space Telescope (HST), 
including a number of circumstellar absorption lines. These ab- 
sorption lines exhibit multi-component structures that are not 
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compatible with a continuous wind, and the derived column den- 
sities indicated a complex differential depletion as a result of 
dust formation. They showed that these structures can be ex- 
plained in principle either by a model consisting of concentric 
shells with different densities, consistent with a highly variable 
mass-loss rate, or by discrete clumps in the line of sight due to 
unknown ejection processes. 

The goal of this work was to construct a more realistic model 
of the a Sco system, including hydrodynamic and plasma effects, 
in order to gain quantitative insight into the dynamics of the cir- 
cumstellar envelope. After a description of our hydrodynamic 
and plasma simulations of the a Sco system and its H n region 
in Sect. [31 we present a comparison to observations in Sect. [3] 



Table 1. System parameters of a Sco used in this work. 



Parameter 


Value 


Source 


T e B,B 


18 200 K 


This work 


loggB 


3.9 (cgs) 


Kudritzki & Reimers (1978) 


M B 


6.7 M Q 


This work 


M A 


18 M Q 


Kudritzki & Reimers (1978) 


d 


185 pc 


Hipparcos catalogue 


6 


23 ±5° 


Reimers et al. (2008) 


Separation (in 2006) 


2'.'73 


Reimers et al. (2008) 


sin ( 


1 


Hopmann(1958) 


e 





Hoomann (1958) 



Notes. S is the position angle of the B star with respect to the primary 
(see text), i and e are the inclination and eccentricity of the orbit, re- 
spectively. 



2. Model calculations 

Our model includes the hydrodynamic effects of the hot H n re- 
gion around the secondary star, a Sco B, on the density and 
velocity distribution of the circumstallar envelope of the binary 
star. The temperature and ionization inside the Hn region are 
calculated by considering both microphysical and radiative pro- 
cesses (see Sect. I2.3l l. and the result is used as input into the 
hydro code (Sect. l2~2l ). 



2.1. Stellar and system parameters 

a Sco is a visual binary system consisting of an Ml lb-type 
red supergiant and a B2.5 V-type dwarf, and i ts distance is d — 
185 pc according to the Hipparcos catalogue. Hopm anr] d 1 95 8l) 
determined the inclination of the orbit and found i ~ 90°. The 
angular separation is 2773 and the B star is ~ 224 AU behind the 
supergiant, which co rresponds to a position angle 5 — 23 ± 5° 
(iReimers et al.ll2008l). 

According to Kudritzki & Reimers d!978l) the effective tem- 
perature of the B star is T e s$ = 18 500 ± 1500 K, and its sur- 
face gravity is loggB = 3.9 + 0.2 (in cgs units). These val- 
ues result from profile fits of Balmer and Hei lines, and the 
mass Mb of the B star was determined by a comparison of 
these values with theoretical evolutionary tracks, yielding Mb = 
7.2 + 0.5 M and Rb _ 5.2 + 1.3 Ry. A pplying more recent evo- 
lutionary models dBressan et al.ll 1 993b gives similar resu l ts, i. e . 
M B = 6.7 ± 0.7 M Q and R B = 4.8 ± 1.1 ^. iBrott et all (l201lh 
presented calculations including rotational effects. Their evolu- 
tionary sequence for 7 M and i> rot = 284 kms~' matches the 
sur face gravity and effective temperature derived for a Sco B 
by [Kudritzki & Reimers d!978l) . who observed a projected ro- 
tational velocity of v mt s in; = 250 kms" 1 . Based o n these re- 
sults and the analysis of iHiellming & Newelll (Qj983), who gave 
a constraint on the Lyman continuum luminosity of the B star, 
we adopted 7/ eff , B = 1 8 200 K. 

For th e mass of a Sco A we adopt ed the value Ma ~ 1 8 M 
found bv lKudritzki & Reimers! d!978l) . Table[T]lists the parame- 
ters used in this work. 



in 3D and uses the equation of state of an ideal gas, 



(r-i)p' 



(i) 



where e is the specific internal energy, p the pressure, p the mass 
density, and y = 5/3. We used the finite-volume i ntegrator that 
is bas ed on a modified Lax-Friedrichs approach (Barmin et al. 
1996), which is implemented in the AMRCART code. This ap- 
proach uses a predictor-corrector procedure to obtain a scheme 
with a second-order accuracy. We did not include the gravita- 
tional potential of the B star, because we were interested in the 
large-scale effects that are caused by the heating that is due to 
its radiation, and the resolution of the computational grid is not 
sufficient for a simultaneous calculation of possible accretion ef- 
fects. The heating by the radiation of the B star is included as a 
source term in the Euler equations (see Sect. I2.31 >. 

2.2.1. Boundary and initial conditions 

We assume that the wind is rapidly accelerated to its terminal 
velocity v x near the supergiant. This assumption does not have 
a significant influence on the results, because the region of in- 
terest, where the wind interacts with the H n region, is far from 
the supergiant, at a distance of the order of ~ 300 AU. Thus, we 
simply apply the terminal velocity at the boundary cells, i. e. 



(2) 



where R& is the radius and »a the orbital velocity of the primary 
star, and d ot points radially outward from the primary. 

At the outer boundaries of the computational domain "free 
flow" is assumed, i. e., density and velocity in boundary grid 
cells equal the values of the cells at the edge of the domain. The 
initial velocity field is given by Eq.[2]at all radii r > R\, and the 
initial density distribution is spherically symmetric with respect 
to the primary, i. e. 

p(r) = 4^T' (3) 

TV I / L/CO 



2.2. Hydrodynamic simulations 

For the hydrodynamic simulations we used the hydro version 
(without magnetic fields) of the AMRCART code , which is 
part of the A-MAZE package dWalder & Folinil l2000) and com- 
prises the adaptive me sh refinement algorithm developed by 
Berg er & Colellal dl989). AMRCART solves the Euler equations 



2.2.2. The computational grid 

As the system is symmetric with respect to the plane of the or- 
bit, we place the center of mass near the upper boundary of the 
computational domain, which reduces the computing time by re- 
ducing the total number of grid cells. The plane of the orbit is 
parallel to the xy plane, and the center of mass of the system is 
placed in the center of that plane. For most of the simulations 
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presented in this work, two levels of refinement are used with a 
refinement ratio of 2 between the respective grid spacings. The 
base grid comprises 160 x 160 x 44 grid cells corresponding to a 
size of a x x a y X a z ~ 12 433 x 12 433 x 3419 AU, where a x , a y , 
and a z are the edge lengths of the domain in x, y, and z direction, 
respectively. 

For the refinement of the grid, the local truncation error 
TE(r, t) of the density is estimated for each cell of the base grid 
with a method that is based on Richardson extrapolation (see 
Berg er & Colellal fT989). If the truncation error of a grid cell ex- 
ceeds a given tolerance limit e, the grid cell is flagged. In the 
next step, cuboid-shaped subgrids are generated that include all 
regions containing flagged cells. We used an error tolerance of 
s = 4 x 10~ 4 and enhanced the resolution by a factor of two in 
the subgrids in our simulations. The error tolerance was applied 
to the relative truncation error, i. e. TE(r, t)/p(r, t) < s. The con- 
tours of the subgrids are indicated as gray rectangles in Fig.Q] 

For the calculation of the heating by the radiation of the B 
star (see next section), spherical coordinates are used, which are 
then converted to the cartesian grid of the hydro code by linear 
interpolation. In Appendix lAl we describe the generation of the 
grid that provides the temperature around the B star as a function 
of position in 3D space. 



S 6 




4 6 8 10 12 
x (1000 AU) 

Fig. 1. Quasi-stationary density and velocity distribution result- 
ing from a simulation with M = 2 x 10~ 6 M yr~' and v x = 
20 kms -1 after ~ 1 orbit or ~ 2560 yrs. The dotted arrow in- 
dicates the line of sight to the B star, whose position is marked 
by the asterisk. The gray rectangles indicate the contours of the 
high-resolution subgrids (see text). 



2.3. Plasma simulations 



The temperature and ionization inside the H n region around the 
B star were calculated wi th version 08 . 00 of the plasma code 
Cloudy, last described by iFerland et alJ (Il998l) . Cloudy uses a 
spherical geometry in ID. In order to obtain a 3D temperature 
distribution of the a Sco system we calculated different models, 
each one corresponding to a different direction, starting from the 
position of the source of ionizing radiation, i. e. the B star (see 
Appendix [A] for details). This approach is justified as Cloudy 
uses the escape probability formalism for the calculation of ra- 
diative transfer effects, so that different points inside the circum- 
stellar envelope do not interact with each other. For the radiation 
of the B star we used a synth etic spectrum given by a n interpo- 
lation on the TLUSTY grid of lLanz & Hubenvl d2007l) . 

The resulting temperature distribution is then passed to the 
hydro code as a source term affecting the total energy. The ra- 
diative heating is assumed to take place at a much smaller time 
scale than the hydrodynamic processes. This procedure is re- 
peated at regular intervals during the calculation, typically every 
other time step as measured on the base grid. The information 
about the state of ionization can be used later in the analysis of 
spectral lines. 

For the treatme nt of Fen, Cloudy provides a large model 
atom developed bv lVerner et al.l dl999i) . It includes the 371 low- 
est levels of the Fe + ion, i.e., up to 93487.65 cm . 

3. Results 

3.1. Density and velocity distribution 

The presence of the hot H n region inside the wind of the su- 
pergiant causes deviations from spherical symmetry. Figure Q] 
shows a cut through the center of mass of the system parallel to 
the plane of the orbit, including the density distribution as well 
as appropriately scaled velocity vectors for a simulation with 
M = 2 x 10~ 6 M yr~' and v x = 20 kms~' after ~ 1 orbit 
or ~ 2560 yrs. At the boundaries of the Hn region around the 
B star that face the supergiant the density increases, and it de- 
creases inside the wake of the H n region. The curved shape of 
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Fig. 2. Enlarged central region of Fig.Q] The contour indicates 
90% ionization of hydrogen. 



the structure results from the orbital motion. The orientation of 
the velocity vectors is still approximately spherically symmetric, 
but the absolute value increases by ~ 50% in the low-density re- 
gions. Figure [2] shows the central region of Fig.Q] including the 
H ii region as a contour that indicates 90% ionization of hydro- 
gen. 

We calculated models with different values of M (0.2, 0.5, 
1, 2, and 5 x 10~ 6 M yr-') and u M (10, 20, and 40 kms" 1 ), 
including all possible combinations of M and v m in these ranges. 

3.2. Comparison to observed circumstellar absorption lines 

Circumstellar absorption lines in t he spectrum of a Sco B have 
been observed wi th GHRS/HST dBaade & Reimersll2007l) and 
with UVES/VLT dReimers et afl 120081) . In Sect. 13.2.11 we de- 
scribe the calculation of synthetic line profiles based on the sim- 
ulations presented above. These synthetic profiles are compared 
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to observations in Sect. l3.2.3l For that comparison, we calibrated 
the wavelength scale of the GHRS spectra by use of the UVES 
spectra (see Sect. I3.2.2I I. It should be noted that the resonance 
lines (0 eV) of abund ant ions are a f fected by interstellar contam- 
ination, as shown by Hagen eTaTI d!987l) . Exceptions are those 
elements which are heavily depleted in the IS medium. 

3.2.1. Synthetic line profiles 

In order to derive the profiles of circumstellar lines we use the 
density and velocity information from the hydrodynamic simu- 
lations to calculate a pure absorption line profile given by 



I v = 7 B exp 



(-£ 



(4) 



where l v is the emergent intensity, 7b the intensity of the B star 
continuum, % v the opacity, and z the coordinate along the line 
of sight towards the B star, which is located at z = Zb- z m ax is 
chosen so large that the line profile does not change if a higher 
value is used for this limit. After about one half of an orbit the 
density structure becomes quasi-stationary and does not change 
its shape anymore. Therefore, the results from the hydrodynamic 
simulations at any time step can be used for the calculation of the 
line profile. We used the density information corresponding to 
the time after about one orbit, i. e. at t ~ P, where P ~ 2560 yrs 
is the orbital period of the system. 

The resulting profiles deviate considerably from the profiles 
resulting from a continuous, spherically symmetric model. As 
an example, Fig. [3] shows the simulated profile of the Zn n UV 
mult. 1 2062.660 A absorption line, based on a simulation with 
M = 10~ 6 M yr~' and u M = 20 km s -1 , along with a plot of the 
same line for the case of an undisturbed, spherically expanding 
wind with the same parameters. Obviously, the absorption com- 
ponent at ~ -20 kms" 1 is present in both lines and coincides 
with the terminal wind velocity, but the absorption profile that 
is based on the hydrodynamic simulation exhibits a more com- 
plex structure, which finds expression in three additional min- 
ima, namely at 29, -7, and 3 km s~' . The radial velocities are 

measured relative to the center of mass of the system. Thermal 
broadening is included by applying a constant temperature of 
5000 K, which is approximately the mean temperature inside the 
H ii region. Microturbulence is not included in the calculation of 
the line profiles, so that all additional effects that influence the 
line profile are of hydrodynamic origin. 

The distinct components seen in the absorption line pre- 
sented in Fig.[3]are related to features of the density and velocity 
distribution along the line of sight to the B star. Figure [4] shows 
the number density of Zn + and the radial velocity as a function 
of distance to the B star. One important result is that the velocity 
is no longer a monotonic function of distance when hydrody- 
namic effects are included, but there is a broad minimum around 
~ 4000 AU. 

In order to test whether the simulations have converged, we 
calculated an additional line profile based on a simulation with 
the same parameters, but with a resolution that is enhanced by 
a factor of two (see Fig. |3j. The result does not show any sig- 
nificant changes. For this high-resolution simulation we used a 
computational domain with a reduced size of a x x a y x a z ~ 
6217 x 6217 x 3419 AU, because it would have been compu- 
tationally too expensive to calculate a high-resolution model of 
the same size as the low-resolution models. This reduced com- 
putational domain does not cover the whole high-velocity region 
that produces the absorption component at 29 kms~'. That 
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Fig. 3. Calculations of the Znn UV mult. 1 2062.660 A absorp- 
tion line based on a hydrodynamic simulation including ion- 
ization with M = 10~ 6 M yr _1 and = 20 kms -1 , and on 
an undisturbed, spherically expanding wind. As a convergence 
test, the hydrodynamic results (low resolution) are compared to 
a model with the same parameters but with a resolution that is 
increased by a factor of two (high resolution). The difference in 
the component at ~ -29 kms 1 is due to the reduced size of 
the computational domain in the high-resolution simulation (see 
text). The abscissa gives the radial velocity relative to the center 
of mass of the system. 
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Fig. 4. Radial velocity relative to the center of mass and Znn 
number density along the line of sight to the B star, based on 
a hydrodynamic simulation and on an undisturbed, spherically 
expanding wind. These data correspond to the absorption line 
profiles shown in Fig.[3](M = 10" 6 M yr _1 , v„ = 20 kms" 1 ). 
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Fig. 5. Absorption profiles of the Znn UV mult. 1 2062.660 A 
and Cm UV mult. 1 2062.236 A lines for different mass-loss 
rates (u^ = 20 km s -1 ). The mass-loss rates are given at the bot- 
tom right of the figure in units of 10~ 6 M yr -1 . 
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Fig. 6. Ionization fractions of different ions as a function of dis- 
tance to the B star along the line of sight. The data are derived 
from a simulation using M = 10 -6 M yr -1 and v m = 20 km s . 

Table 2. Velocity shifts applied to the GHRS spectra (see Fig . [7J . 



Transition 


Shift (kms- 1 ) 


CuiiUVmult. 3 1358.773 A 


-3.2 


Nin 1393.324 A 


-4.2 


Znn UV mult. 1 2062.660 A 


-3.0 


Cm UV mult. 1 2062.236 A 


-3.0 


Fen UV mult. 1 2622.452 A 


-1.6 



means that the broad minimum in the radial velocity (Fig. |4j is 
not fully covered, which results in a weaker absorption compo- 
nent. 

Figure [5] shows the dependence of the Zn n UV mult. 1 
2062.660 A and Cm UV mult. 1 2062.236 A lines on the mass- 
loss rate. While the position of the component at 20 km s -1 

remains nearly unchanged, the position of the other strong com- 
ponent, which is associated to the increased density at the bound- 
ary of the H n region, shifts to lower velocities as the mass-loss 
rate decreases and the H n region becomes more extended. There 

is also a component at 29 kms -1 , which is very weak in the 

Zn ii line. The profiles of the Cr n line show that this component 
is considerably stronger at lower mass-loss rates. This is because 
the wake of the H n region is more extended if the overall den- 
sity is lower, so that the minimum of the radial velocity, shown 
in Fig. |4]for M = 10 -6 M yr -1 , becomes broader and a higher 
number of particles along the line of sight move with veloci- 
ties larger than the terminal wind velocity o„o. That is, the larger 
extent of the volume containing gas at high velocities, which 
produces the absorption component at -29 km s , overcompen- 
sates the effect of the decreased density. 

In the Cr n line there is practically no absorption at positive 
velocities, because most chromium is doubly ionized throughout 
the H ii region. In contrast, Zn ++ is only present close to the B 
star, which is partly due to the different ionization potentials of 
Zn + (144893 cm' 1 ) and Cr + (132966 cirT 1 ). The relative abun- 
dances of H + , Cr + , and Zn + are shown in Fig. [6] as a function of 
distance to the B star along the line of sight. 



3.2.2. Calibration of the observed wavelength scale 

The positions of the observed absorption components in 
the GHRS spectra ar e slightly different in every line 
(Baade & Reimers 2007), which is due to the relatively large er- 
ror of the wavelength calibration. The data were recorded using 
the Large Science Aperture of GHRS, which results in a maxi- 
mum error of ~ 4.5 km s _1 dHeap et al]|1995l) . The wavelength 
calibration of the spectra observed with UVES is much mo re ac- 
curate. It amounts to > 120 ms -1 dD'Odorico et al.ll2Q00t) . The 
panel at the top right of Fig. |7]shows the Ca n mult. 1 (H and K) 
lines. There are three components in both lines, at -21.0, -13.6, 
and -5.9 kms -1 , respectively. 

Now, supposing that the component at -21 kms -1 in the 
Can lines is associated to the terminal wind velocity as sug- 
gested by the synthetic line profiles presented in the previous 
section, the positions of the GHRS lines can be calibrated by 
applying shifts that place the most blue-shifted strong compo- 
nent at -21 kms -1 , which is acceptable as long as the shifts are 
smaller than the maximum error of the GHRS wavelength cal- 
ibration (4.5 kms"" 1 ). Table [2] lists the shifts corresponding to 
the five GHRS lines presented in Fig. [7] The Cm line is satu- 
rated, so that the positions of the strong absorption components 
are rather uncertain. Therefore, the shift determined for the Zn n 
line is used, which lies very close to the Cm line. For the Fen 
line the central component is used to define the shift, because the 
observed component near -21 kms -1 is weak and not clearly 
pronounced. 

Observed absorption lines of Ti n also show a strong absorp- 
tion component at ~ -21 kms -1 , which may serve as an addi- 
tional justification for identifying the position of this component 
with the terminal wind velocity. The two panels at the bottom 
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Fig. 8. Comparison of theoretical line profiles (solid curves) cor- 
responding to a mass-loss rate of M = 2 x 10~ 6 M yr _1 and 
a wind velocity of Voo = 20 km s with observations (dashed 
curves), including a velocity calibration provided by the UVES 
spectra (see Sect. 13.2.21 ) and a systemic velocity of - 1 .3 km s . 
The relative flux is given as a function of the observed radial 
velocity. 

right of Figure Q show absorption line profiles of the Ti n multi- 
plets 1 and 2. 

3.2.3. Observed line profiles 

Figure [7] shows a selection of absorption lines observed with 
GHRS/HST and UVES/VLT in the spectrum of a Sco B. The 
shapes of the line profiles produced by different ions or different 
energy levels of the same ion (Tin) show a large variety. The 
observed absorption features depend on the excitation of fine- 
structure levels of the ions and the ionization of the correspond- 
ing element along the line of sight. Most of the lines exhibit three 
strong absorption components at the approximate positions -21, 
-14, and -6 kms , respectively. Some of the lines produced 
by the eV levels have an additional component at -29 km s 
(Tin 3383.759 A, Cm 2062.236 A, Znn 2062.660 A, and Nin 
1393.324 A), which can be explained by the broad minimum of 
the radial velocity in the wake of the H n region far from the B 
star (see Figs.|5]and|4]l. 

For the comparison presented in the following, a systemic 
velocity of u sys = —1.3 kms -1 is added to the wavelength scale 
of the simulated profiles, which results from a comparison of 
the projected orbital velocities resulting from the adopted orbital 
confi guration (6 — 23°, e — 0) with the observed radial velocities 
(see lReimers et al.ll2008l) . 

Cm, Cu ii, Nin, and Znn Figure [8] shows a comparison of the 
observed Cm, Cun, Nin, and Znn lines presented in Fig. [7] to 
theoretical profiles corresponding to a mass-loss rate of M — 
2xl0~ 6 M Q yr~' and a wind velocity of v x - 20 km s _1 . The sim- 
ulated line profiles are convolved with a gaussian profile function 
with a FWHM of 3.5 km s _1 to account for the limited resolution 
of the GHRS/HST spectra. 

Obviously, there are a number of differences between the 
simulated and the observed line profiles. First of all, there are 



only two strong components in the simulated profiles, whereas 
there are three in the observed profiles. Secondly, the observed 
component at -6 km s is much stronger than predicted by the 
model calculations in the Znn and Cun lines. However, with 
the wavelength calibration described in Sect. I3.2.2l the positions 
of the two strong components in the simulated profiles approx- 
imately match the positions of the outer two of the observed 
strong components. It turns out that, within the error of the wave- 
length calibration, the observed GHRS profiles are consistent 
with the assumption that the position of the most blue-shifted of 
the strong absorption components is at -21 kms -1 , as suggested 
by the absorption lines measured with UVES and by the theoret- 
ical lines, all of which include a strong absorption component at 
the position corresponding to the terminal wind velocity. 

The weaker absorption component at ~ -29 kms 1 , which 
is clearly visible in the Cr n and the Zn n line, is very weak in 
the corresponding theoretical profiles. The simulations predict a 
stronger absorption component at this position for lower mass- 
loss rates, especially in the Cm line (see Fig . |5j ■ This may be an 
indication of a time-dependent mass-loss rate, which could pro- 
duce distinct absorption components corresponding to different 
mass-loss rates (see Sect. I4.21 i. 

The equivalent widths of the observed profiles and the simu- 
lated profiles shown in Fig.[8]agree well. Moreover, the compo- 
nent at 21 kms 1 of the simulated profiles gives a good fit 

of the observed components. The absorption in the components 
at -6 km s is stronger in the observed lines, which may be due 
to interstellar absorption as all the lines result from transitions 
starting from the lowest energy levels of the ions. The theoreti- 
cal absorption lines resulting from simulatio ns using a mass-loss 
rate of M = 10~ 6 M yr _1 , as given by iReimers et al.l (2008) 
based on a simplified model of the circumstellar shell neglecting 
dynamic effects, are much weaker than the observed ones. That 
shows the significant impact of the hydrodynamic effects, which 
will have to be considered in future studies of mass loss. 

Differences between the observed and simulated profiles are 
to be expected, because the simulations do not account for a 
possible time dependence of the mass-loss rate. Indications of 
multiple shell ejection have been observe d in IR images in both 
a Sco and a Ori (e. g. iDanchi et al.ll 19941 see also Sect.Q}. The 
additional third component in the observed lines seen in Fig.|8]at 

14 km s could well be the result of such a discrete ejection. 

An additional complication of a comparison o f observations with 
theory is the observed differential depletion dBaade & Re imers 
2007). These effects may be the reason for a more complicated 
structure of the absorption profiles. 

Fe ii The observed absorption profile of the Fe n UV mult. 1 
2622.452 A line presented in Fig. Qresults from the fine struc- 
ture level in the ground term at 977.053 cm -1 . The population 
of this level cannot be reproduced by applying a constant scale 
factor to the total Fe + density, as would be expected if the levels 
were populated according to their statistical weights. 

The model atom of IVerner et al.l (1 19991) that was used in 
the plasma simulations yields the population of the level at 
977.053 cm -1 , which is the highest level in the ground term. 
It turns out that the relative population of the level is very low 
inside the cool neutral wind. The population of the Fe n energy 
levels depends on the electron temperature, which determines 
the rate of collisions. The temperature inside the Hn region is 
~ 5000 K, but the temperature in the Hi region is not exactly 
known. 
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Fig. 7. Absorption profiles observed in the spectrum of a Sco B with HST/GHRS (Cu n, Fe n, Ni n, Zn n, and Cr n lines) and 
VLT/UVES (Ti ii and Ca n lines). 



Figure [9] shows the absorption line profiles corresponding to 
two different lower limits of the electron temperature along with 
the observed absorption profile, which is shifted by -1.6 km s 
according to Table [2] (cf . Sect. I3.2.2l i. In these simulations, the 
temperature was not allowed to fall below a given lower limit. 
Obviously, the relative population of the 977.053 cirT 1 level 
steeply decreases after a peak at the boundary of the He re- 
gion. The slope of the right flank of the strong absorption com- 
ponents in the theoretical profiles does not exactly match the 
observed one, because there is a considerable amount of reemis- 
sion, which is not included in the calculation of the theoretical 
line profiles. The decrease of the absorption towards lower ve- 
locities qualitatively matches the observation, although the slope 
is much too steep in the theoretical profiles. The absence of 
the middle absorption component in the theoretical profiles is 
not surprising since it is absent in all theoretical profiles (cf. 
SectEZLQ. 

The assumption that the neutral part of the wind has a tem- 
perature of 4000 K is probably not realistic, and even in this case 
the predicted absorption is too weak in the outer parts of the en- 
velope. A value of T % 1000 K in the H i region is more realistic. 
Apparently the calculation of the level population with Cloudy 
does not give correct results, maybe because of an underestima- 
tion of the rates of collisional or radiative processes that populate 
the level that is responsible for the absorption, or because of ad- 
vection effects, which are not included in the simulations (see 
Sect. |4|. In the plasma simulations, the gas is assumed to be at 
rest with respect to the B star, which may introduce another un- 
certainty in the calculation of the populations of fine-structure 
levels, as the velocity field may have considerable effects on 
the continuum pumping rates. However, the strength of the ab- 
sorption component at 6 kms -1 approximately matches the 

strength of the observed component, so that the mass-loss rate of 
2 X 10" 6 M Q yr _1 , which is consistent with the observed profiles 
of the lines of other ions (cf. Fig. [8), can be confirmed. 

The results presented in this section suggest that the weak- 
ness of the observed absorption components of this Fe n line at 
-14 and -21 kms~' is due to the relatively high energy differ- 
ence between the lowest level and the fine-structure level corre- 
sponding to the transition. Apparently the mechanisms that ex- 
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Fig. 9. Observed and theoretical absorption profiles of the Fe n 
UV mult. 1 2622.452 A line for M = 2 x 10~ 6 M yr' and 
Voa = 20 km s . For the theoretical profiles, two different lower 
limits of the electron temperature were used in the correspond- 
ing plasma simulations (see text). The relative flux is given as a 
function of the observed radial velocity. 



cite Fe n to higher fine-structure levels are not efficient in the 
outer parts of the circumstellar envelope, i. e. at low densities, 
far from the B star and its H n region. 

Ca ii and Ti n The number densities of Ca + and Ti + along the 
line of sight cannot be reproduced in the framework of the 
present model. In the case of Ti + , the plasma simulations yield 
a density that is much too high and density gradients that are 
too steep to reproduce the observed absorption profiles. The ion- 
ization of Ti + is complicated because its ionization potential 
(13.576 eV) is close to the Lyman edge. Simulations show that 
the ionization of Ti + sensitively depends on the local temperature 
in the circumstellar shell. The absence of the absorption compo- 
nents at -6 and -14 km s in the observed line profiles of Ti n 
shown in Fig. [7] suggests that most titanium is at least doubly 



7 



K. Braun et al.: A hydrodynamic study of the circumstellar envelope of a Scorpii 



ionized up to a large distance to the B star, i. e. even far outside 
the H ii region. However, the absence of these components may 
also be due to dust depletion. The component at -30 kms~' is 
observed only in the eV transitions of Ti n, which supports evi- 
dence from the simulations that the line is formed far away from 
the system. In the case of Ca + , the predicted number density is 
much too low. In the whole circumstellar envelope Ca ++ appears 
to be the dominant stage of ionization of calcium, and Ca + shows 
a very complex dependence on the varying conditions along the 
line of sight. Thus, the Ca n lines are inappropriate for mass-loss 
diagnostics. 

The simplified treatment of the corresponding model atoms 
may be the reason for the discrepancy between theory and obser- 
vation in the case of Ca n and Ti n. The Cloudy code treats Ca + as 
a five-level atom, which includes the fine-structure levels of only 
the lowest three terms (4s, 3d, and 4p) with the corresponding 
transitions (multiplets 1, 2, and IF). Ti + is treated in the frame- 
work of the scandium-like isoelectronic sequence as a five-level 
atom without resolving any fine-structure components, includ- 
ing only the terms a 4 F, z 4 G°, z 4 F°, and the terms 4 F° and 4 D° 
of the configuration 3d( 2 D)4s4p( 3 P°). Thus, for Ti + , populations 
of individual fine-structure levels are not calculated. Moreover, 
only transitions to the ground term of Ti n are included in the 
simulation. The depletion of titanium due to dust formation is 
probably the most signi ficant effect produci ng deviations from 
the observed data (cf. Baade & Reimers 2007). 



3.3. Comparison to optical line emission from the Antares 
nebula 

The Antares emission nebula, which is associated to the H n re- 
gion around a Sco B, is seen in the optical in Ha and various 
other emission lines, last described in detail by Reimer s et al.l 
(2008). The simulations presented above will be compared to 
observations of the Ha emission in Sect. 13.3.11 Moreover, for- 
bidden transitions of Fe n, which produce the most prominent 
emission lines of the Antares nebula, will be used to compare 
observations to our model in Sect. 13.3.21 The observations were 
performed with UVES/VLT using a long slit perpendicular to 
the line connecting the two stars, and yielded a scan of the emis- 
sion nebula with a step size of ~ O'.'S, starting 079 west of the 
supergiant (see Fig.fTOb. 



3.3.1 . Ha emission 

The plasma simulations show that the strongest Ha emission 
in the Antares nebula emerges near the B star, where most hy- 
drogen is ionized. Outside the H n region the emission is much 
weaker and follows the total density (Fig. [TJ. In an Hn region, 
Ha is naturally produced by recombination of electrons and pro- 
tons. Another important mechanism is pumping by Lyman line 
photons from the source of radiation, so that the presence of 
Lyman absorption or emission lines in the spectrum of the source 
of radiation affects the production of Ha phot ons. The B-star 
atmos pheres presented in the TLUSTY grid (lLanz & Hubenvl 
|2007|) show strong absorption in the Lyman lines. As the Antares 
H ii region is optically thin in Ha, the observed emission is solely 
determined by the population of the upper levels and the oscil- 
lator strengths. The sensitivity of the Ha emission produced in 
the H ii region to changes of the input flux in the range of the 
Lyman lines requires a high-resolution input spectrum in the 
Cloudy simulations. For the calculations presented in this sec- 
tion we used a resolution of Av/v = 5 X 10~ 4 . 




Fig. 10. Slit positions used for the observations of the Antares 
nebula. The rectangles indicate the position and size of the slit 
in the blue arm relative to the a Sco s ystem as it is seen in the 
sky. This is Fig. 1 from Reimers et alj d2008l) . 



The spatial extent of the Hn region is determined by the 
mass-loss rate of the supergiant and the Lyman-continuum flux 
of the B star, and can be measured by observing th e spatial extent 
of the Ha emission from the Antares nebula (cf. iReimers et al.l 
2008). Figure[TT]shows the Ha intensity distribution at 3'.' 4 from 
the supergiant, measured on the projected line connecting the 
two stars (see Fig. [Tol l. The shape and extent of the observed 
and theoretical Ha intensity distributions agree, but the center 
of the theoretical distribution is shifted to higher velocities with 
respect to the observed distribution. This is probably due to the 
uncertainties in the geometry of the system and can be explained 
by an overestimation of the position angle of the line connecting 
the two stars relative to the plane of the sky (e = is only a 
rough estimate, and due to the long period it will be difficult to 
improve on this). 

The intensity in Fig. QT| is given as a function of velocity. 
Assuming pure emission, the emergent intensity is given by the 
integral 



IM = j s ' 



j(s)®[v r (s) - v]ds, 



(5) 



where s ranges from s m ; n = -0.3a to s max = +0.3a along the line 
of sight, and the zero point lies in the plane defined by the line 
connecting the two stars and the direction perpendicular to the 
orbit, so that the H n region is fully included, a is the edge length 
of the domain in x and y direction. s m j n and s max were chosen 
such that the integration includes a maximum of the available 
simulated data but does not exceed the boundaries of the com- 
putational domain at the outermost slit position. As most of the 
Ha emission comes from the H n region, this range ensures that 
the integration covers the whole emitting region, j is the emis- 
sion coefficient for isotropic emission, 



hv lu n u A lu 
An 



f 

Jo 



(6) 



resulting from the frequency vi u of the transition, the number 
density n u of ions in the upper state, the transition probability 
Aiu, and the profile function <p(v). h is Planck's constant, v r (s) 
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integration of the theoretical data reads 
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Fig. 11. Ha emission at 3'.' 4 from the supergiant. The top panel 
shows the observed flux in arbitrary units, the bottom panel the 
intensity derived from a simulation with M — 2 x 1CT 6 M Q yr _1 
and Voc = 20 kms -1 . The ordinate indicates the position along 
the slit. As can be seen, cf. also Fig. [12] the theorical distribution 
is slightly too narrow, which means that the adopted mass-loss 
rate is here slightly too large. 



the radial velocity at s, and O a gaussian profile that introduces 
thermal broadening, 



<S>(v) = 



1 



Add y/n 



exp 



(Au D ) 2 



(7) 



where the Doppler width Auu is defined by pure thermal broad- 
ening with a temperature of 5000 K. Microturbulent broadening 
is not included, so that the hydrodynamic effects can be clearly 
identified. /„ is the intensity per unit radial velocity, and the in- 
tensity per unit frequency reads I v = AqI u , where Ao is the rest 
wavelength of the transition. 

For a comparison with the observations as presented in 
Fig. QT| for Ha, the spectral resolution of UVES and the see- 
ing have to be considered. An analysis of the lines of the wave- 
length calibration lamp yields an average FWHM of 3.5 kms 1 
in the vicinity of the [Fen] mult. 20F 48 14.55 A line, which we 
present in the next section, in the range from 4789 to 4832 A. We 
adopt this value for the spectral resolution of the UVES spectra. 
The seeing during the observations was ~ 076. Therefore, the 
theoretical intensity distributions are convolved with a gaussian 
profile of FWHM = 3.5 kms 1 along the frequency coordinate 
and with another gaussian corresponding to the seeing in the di- 
rection of the spatial coordinate along the slit. A systemic veloc- 
ity of -1.3 kms~' is added to the theoretical velocity scale (see 
Sect. 1372731 . 

For a comparison of the spatial extent of the observed and 
theoretical emission, it is easier to compare the frequency- 
integrated emission as a function of position along the slit. The 



£ 



j(s)ds. 



(8) 



The results are shown for different slit positions in Fig. [T2lfor 
a simulation with M = 2 x 10~ 6 M yr _1 and = 20 kms~', 
along with the corresponding observational data. In this simula- 
tion, the slit skims only the outermost part of the H n region at the 
slit positions 179 and 5'.'4. Therefore, the theoretical data at these 
slit positions should be interpreted with care, because outside the 
H ii region the Ha production is dominated by line pumping ef- 
fects, which may not be well reproduced with the simplified ap- 
proach used by Cloudy for the radiative transfer (see Sect. 12. 31 . 
The central minimum in the observations corresponding to 274 
is probably due to an artifact related to the data reduction. At 
the other slit positions, the simulated and observed distributions 
agree well in shape and extent. 

At 274 and 279, the theoretical Ha profiles are narrower than 
the observed profiles, indicating a lower mass-loss rate, while 
the observed data at the other slit positions further away from 
the supergiant are better reproduced by the model. This may be 
due to the mass loss rate being time-dependent. 

3.3.2. [Fe n] emission 

The most prominent emission lines seen in the Antares neb- 
ula are the forbidden iron lines. The spatial and spectral distri- 
bution of th e fine at 4814.55 A ( multiplet F20) was presented 
in detail by Reime rs et a D (120081) . The data obtained with the 
Cloudy calculation s, which include the large Fe + model atom of 
IVerner et alj (H999D (see Sect. 12.31 ). can be used to compare the 
results of the simulations to the observed data. 

As an example, Fig. Q~3] shows the distribution of the [Fen] 
mult. 20F 4814.55 A intensity at 374 both for the observed and 
the simulated data. Obviously, the simulations do not yield a re- 
alistic picture of the [Fe n] emission for that slit position. The 
simulations suggest a circular structure around a central maxi- 
mum at ~ 5 kms -1 , while the maximum flux in the observed 
data is concentrated in a more compact structure, which is ap- 
proximately circular but open to the bottom, between ~ and 
~ 8 kms 1 as measured at the center of the slit. The overall 
extent of the observed emission approximately agrees with the 
theoretical distribution and is consistent with the extent of the 
Ha emission in both the velocity and the spatial coordinate. 

The discrepancies between the observed and theoretical 
[Fe n] emission could be due to an inadequate treatment of the 
Fe + ion in the Cloudy code. The emission may also be influ- 
enced by advection effects, which may be significant especially 
at the western boundary of the H n region towards the open side 
of the wake, where the density is small. An alternative expla- 
nation would be that due to dust depletion in the high density 
boundaries of the structure, a large fraction of Fe is not available 
for Fe ii emission. This effect seems to increase with distance to 
the B star. 

A common feature of the observed forbidden and allowed 
Fe + emission is an asymmetry with respect to the center of the 
slit, which is clearly visible e. g. at 574 from the supergiant 
(Fig. n~4b . Apparently, there is more emission from the north- 
ern than from the southern half (see also Fig. [13). This is an- 
other indication of density structures that are not included in the 
model, which is symmetric with respect to the plane of the orbit. 
A deviation of the inclination from 90° would also cause such an 
asymmetry. However, the Ha flux distributions do not show any 
significant asymmetries. 
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Fig. 12. Frequency-integrated Ha flux as a function of position along the slit for a simulation with M = 2 x 10 M yr and 
= 20 km s (solid lines), and the corresponding observations (dashed lines). 
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Fig. 13. [Fen] mult. 20F 4814.55 A emission at 3'.'4 from the 
supergiant. The top panel shows the observed flux in arbitrary 
units, the bottom panel the intensity derived from a simulation 
with M = 2 x 10~ 6 M yr _1 and v M = 20 kms -1 . The ordinate 
indicates the position along the slit. 



4. Discussion 

4.1. The mass-loss rate 

The hydrodynamic simulations presented in Sect. l3.1l show that 
the hot H n region that is moving with the B star through the cir- 
cumstellar envelope of a Sco produces strong deviations from 
spherical symmetry in the density and velocity distributions. The 
theoretical absorption line profiles derived from these simula- 



4 
2 

■I 

c 
a, 

-2 
-4 



■I 

6 
a, 

-2 
-4 





P 



0.2 



0.1 



0.0 



1.0 




0.8 




X 






0.6 






1 


0.4 


pit 






2.0 


o 






0.0 





-30 -20 -10 10 

Radial velocity (km s _1 ) 



20 



Fig. 14. [Fen] mult. 20F 4814.55 A emission at 5'.'4 from the 
supergiant. The top panel shows the observed flux in arbitrary 
units, the bottom panel the intensity derived from a simulation 
with M = 2 x 10~ 6 M Q yr _1 and = 20 kms. The ordinate 
indicates the position along the slit. Obviously, there is more 
emission from the upper (northern) half in the observed flux dis- 
tribution. 



tions (Sect. 13.2. U exhibit a pronounced multi -component struc- 
ture comparable to the observed line profiles presented in Fig. [7] 
The best match between theoretical and observed line profiles is 
achieved with a mass-loss rate of M = 2 x 10~ 6 M yr _1 and 
a wind velocity of v x = 20 kms -1 (Sect. 13.2.3b . This value 
of the mass-loss rat e is twice as high as the rate derived by 
Reim ers et al.l (120081) . which was based on a spherically expand- 
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ing model of the circumstellar shell. This discrepancy is due to 
the decreased density in the wake of the Hn region (see Figs.Q] 
and|2]i. 

The range of values of M and Uoa that is covered by the sim- 
ulations is rather small, owing to the limitations imposed by the 
large amount of computing time that the hydro code requires. 
As the AMRCART code is not adapted to modern parallel com- 
puting clusters, the hydrodynamic simulations were carried out 
on a single CPU core, while only the photoionization/radiative- 
transfer calculations were executed in parallel. A finer grid in M 
and Voo would improve the accuracy of the resulting mass-loss 
rate, and it might be instructive to calculate models with differ- 
ent semi-major axes and position angles to account for the un- 
certainties in the orbital parameters (see next section). However, 
the theoretical data presented in this work and the available ob- 
servations agree fairly well, and the remaining systematic errors 
(see below) are probably larger than the error that is due to the 
coarseness of the grid of models, i. e. the relatively large step size 
in M and v x . We presume that the total error of the mass-loss rate 
including systematic uncertainties is ~ 5 X 10~ 7 M Q yr _1 . 

An important systemati c error is introduced by the observed 
differential dust depletion dBaade & Reimersll2007l) . This effect 
depends on the considered element and on the local conditions in 
the circumstellar shell. It results in the observed absorption pro- 
files being weaker than predicted for a given mass-loss rate. In 
contrast, interstellar absorption can make the absorption com- 
ponents near kms 1 stronger than predicted by the model, 
which is probably responsible for the discrepancies in the Zn n 
2062.660 A and C u ii 1358.773 A lin es (Fig. [Hi. According to 
ISnow et all d!987l) and lCardellil dl984l) . dust depletion is of mi- 
nor importance for zinc, which means that the mass-loss rate can 
best be determined by use of the Zn n absorption lines. 

In the present approach, the radiative transfer is treated in 
a simplified manner with the escape probability formalism (cf. 
Sect. 12.3b . This is probably the major source of error in the cal- 
culation of the [Fe n] line emission, and the exact solution of the 
scattering problem would probably resolve some of the apparent 
discrepancies between the theoretical predictions and the obser- 
vations. Moreover, some of the circumstellar absorption lines in 
the spectrum of a Sco B, e.g. Cm 2062.236 A, have P Cyg- 
type profiles, which cannot be reproduced with a pure absorption 
model. For a calculation of the reemission that is superposed in 
these lines on the absorption profile, exact radiative transfer sim- 
ulations have to be included in the model. 

4.2. Asymmetries and time-dependent effects 

The observations indicate an asymmetric, non-stationary charac- 
ter of the circumstellar shell of a Sco. The observed absorption 
lines exhibit an additional component at - 14 km s , which can- 
not be reproduced with the assumption of a constant mass-loss 
rate. Moreover, the integrated Ha flux (Fig. [TZl i also appears to 
indicate slightly different mass-loss rates at different distances to 
the B star (Sect. [3~3~TT >. The observations of [Fen] emission re- 
veal that the density distribution is not symmetric with respect to 
the plane of the orbit, which cannot be reproduced in the frame- 
work of the present model. 

The uncertainties in the geometry of the a Sco system may 
partly account for the asymmetry of the observed [Fe n] emis- 
sion and the shift of the theoretical emission in Ha and [Fe n] 
4814.55 A. The orbital parameters, such as the inclination, ec- 
centricity, and the orbital velocities, are not well known because 
of the large orbital period (~ 2560 yrs). Especially the radial ve- 



locity of the supergiant is hard to measure, because the observed 
velocity results from a s uperposition of pulsations and orbital 
motion (ISmith et al.lll989l) . 

The plasma calculations used for determining the temper- 
ature and ionization structure of the circumstellar envelope are 
based on the assumption that the time scales of cooling and heat- 
ing are small compared to the dynamic time-scale. This is not ex- 
act at the western boundary of the H n region, towards the low- 
density wake, where also advection effects may have a signifi- 
cant impact on the ionization balance. A model including time- 
dependent plasma-effects may considerably improve the under- 
standing of the observed absorption and emission features and 
could be the subject of future studies. 

5. Conclusions and Outlook 

The results of the combination of observations with hydrody- 
namic and plasma simulations presented in this work show that 
the circumstellar envelope of the a Sco system is strongly influ- 
enced by dynamic effects. A calculation of absorption line pro- 
files based on the simulated density and velocity distributions 
and ionization structure, and a comparison to GHRS/HST and 
UVES/VLT spectra reveal that the mass-loss rate was underesti- 
mated by a factor of two in earlier studies that were based on an 
undisturbed, spherically symmetric circumstellar shell. The re- 
sulting mass-loss rate of M = 2 x 10~ 6 M G yr~' is confirmed by 
an analysis of the Ha emission from the Antares nebula, which is 
based on spatially resolved emission distributions observed with 
UVES. 

The theoretical absorption line profiles exhibit a multi- 
component structure as a natural result of the influence of the 
hot H ii region that is moving with the B star through the wind of 
the primary. Three of the four observed components can be ex- 
plained accordingly. However, the origin of the additional com- 
ponent at -14 kms~' seen in most absorption lines (Pig. O re- 
mains uncertain and is probably an indication of time-dependent 
mass-loss. 

The structure of the [Fe n] line emission of the Antares neb- 
ula as observed with UVES cannot be reproduced correctly. A 
more sophisticated model including exact radiative transfer cal- 
culations and time-dependent simulations of the ionization bal- 
ance allowing for advection effects would help to understand 
the remaining discrepancies between theory and observations. In 
this context, it would be desirable to parallelize the hydro code, 
which will make it possible to calculate a more extensive grid 
of models. Future projects may deal with these improvements of 
the model calculations. 

A number of empirical mass-loss rates of f Aur systems were 
determined on the basis of sp herically symmetric density and ve - 
locity distributions (see e. g. lChe et alJl983HBaade et al.lll996l) . 
Some of these systems, e. g. f Aur and 31 Cyg, also contain Hn 
regions, which may lead to even more severe dynamic effects 
than in a Sco due to the much smaller orbital periods. The cur- 
rent values of the mass-loss rates of these stars may have to be 
revised on the basis of more realistic models including dynamic 
processes as presented in this work for a Sco. 
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As the Cloudy simulations are independent of each other, 
they can be safely executed on different CPUs. The simulations 
presented in this work are symmetric with respect to the orbital 
plane so that it is sufficient to cover only the lower half of the 
domain, i. e. 

N-2 
i=(JV-l)/2 

For most of the simulations presented in this work, we used 
N = 9,i. e. N^T" = 45. The only exception is the high-resolution 
simulation (cf. Fig. [3}, where TV = 13 and AT™ 1 ™ = 101. 



Appendix A: The temperature grid 

For the calculation of the 3D temperature distribution, the ID 
models calculated with Cloudy are distributed homogeneously 
in 3D space, i. e. with a constant increment in 9. The B star is the 
starting point for all ID models, which give the physical con- 
ditions along a straight line in radial direction. The angle is 
measured from the upper pole, i. e. from the positive z direction 
in the AMRCART coordinate system, downwards, and the an- 
gle <p is measured counterclockwise from the direction towards 
the supergiant. The increment in <p in the plane of the orbit, i. e. 
at 9 = 7r/2, equals the increment in 9, which is determined by a 
given number TV according to A9 = n/(N- 1). Only odd numbers 
are used for TV so that there are always two models describing the 
temperature distribution on the line connecting the two stars, one 
away from the primary star, one towards it. In directions nearer 
to the poles, A<p is increased according to A(f> = 7r/(Nsm9 - 1), 
where the expression TV sin 9 is rounded to the nearest integer that 
is > 2. For f? = and 9 — n one direction (<p = 0) is calculated. 

Thus, the total number of Cloudy simulations computed to 
cover the whole domain for a given N is 

N-2 

N a = 2 + ^2(^-1), (A.l) 

!=1 

where 2(W,- - 1 ) is the number of directions calculated for a given 
9i = iA9, and 

ftfrinft, sin ,.>2//V_ 
12, else 
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